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Résumé 

Nous élaborons un modèle de propagation d'ondes électromagnétiques dans un mi- 
lieu, inliomogène, avec une couche rugueuse, et dissimulant un objet. Nous obtenons 
un milieu effectif, puis nous résolvons le problème par équations intégrales. Pour citer 
cet article : J.-B. Bellet, G. Berginc, C R. Acad. Sci. Paris, Ser. * *** (2010). 

Abstract 

Electromagnetic model of a hidden object. We model propagation of electro- 
magnetic waves in a médium, which is inhomogeneous witli a rough layer, and wliicli 
hides an object. We first get an effective médium, and then we solve tlie problem 
by intégral équations. To cite this article: J.-B. Bellet, G. Berginc, C. R. Acad. Sci. 
Paris, Ser. * *** (2010). 



Introduction 



Nous sommes motivés par l'imagerie laser d'un objet dans un milieu fluc- 
tuant. Ceci est adapté pour la détection d'une tumeur enfouie dans le derme, 
ou d'un engin explosif improvisé camouflé. Nous résolvons un problème mo- 
dèle de propagation des ondes électromagnétiques dans un milieu contenant 
un demi-espace avec des inhomogénéités, surmonté d'une couche rugueuse. La 
couche peut représenter la couche cornée de la peau, ou un fllet de camouflage. 
Le demi-espace peut constituer l'intérieur de la peau, ou un milieu de camou- 
flage. Nous nous plaçons dans un contexte bi-dimensionnel périodique. Enfln, 
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nous considérons un cas asymptotique multi-échelles : la couche est mince par 
rapport à la longeur d'onde, et les inhomogénéités sont petites par rapport à 
la couche mince. 

Ainsi, on modélise la couche supérieure par la couche mince périodique Qi'^ = 
{(xi,X2) : < a;2 < ^f{xi/^)}, où / > est une fonction 1-périodique, et 
> est un petit paramètre (par rapport à la longueur d'onde). Cette couche 
mince est délimitée par les interfaces 70 = {x2 = 0} et 7ç = {{xi,^f{xi/^)),Xi G 
M}. La sous-couche est le demi-espace ^~ = {x2 < 0}. Le milieu extérieur 
est = {^2 > (3^1/0}- Le milieu supérieur et la couche sont 
homogènes. Le milieu inférieur est un milieu ^F-périodique contenant des 
petites inclusions, où la cellule de référence Y = (0,^i) x (0,£2) {h, h > 0) 
contient une inhomogénéité B. Bien que la taille caractéristique de la période 
soit encore notée ^ dans nous la supposons être à une échelle inférieure à 
la taille de la couche. Dans la cellule Y, la perméabilité et la permittivité sont 
respectivement : /iy(|/) = + /i]Iy^;g(î/), eyiy) = -F elyy;^^/), 

où fiB, f^,£B,£ > sont des constantes. Dans tout le domaine M^, la perméabi- 
lité et la permittivité sont respectivement : = yU"'"l^+(x) + /i^^l^ci(x) + 

/iy(x/,^)l^-(x), e^{x) = e~^l^+{x)+e''^lQ,ci{x)+eY{x/^)l^-{x), où fi~^,e~^ > 
sont des constantes et où l'on a prolongé par périodicité les coefficients dans 
^~ , et 1. est la fonction indicatrice. On illumine dans la partie supérieure 
du domaine par un laser, i.e. une onde plane Ui^dx) = e^*^^^'^. Cette onde est 
de pulsation a; > 0, de nombre d'onde = a;A/ë+/ï+, d'angle d'incidence 



e = (6'i,6'2) avec 6*2 < et 



= L La longueur d'onde est := 27r/A;+, 
supposée grande par rapport à la taille de la couche mince ^. Le problème 
modèle consiste alors en l'étude de la propagation de l'onde résultante mç dans 
le milieu : 

V ■ — Vuç + w^eçuç = 0, dans (1) 

avec en plus une condition de radiation. Cette dernière sélectionne la solution 
du physique du problème ; en particulier, elle veille à ce que l'onde réfléchie 
— Mine dans le milieu supérieur et l'onde transmise dans le milieu 
inférieur soient des ondes sortantes. Enfin, une tumeur dans le derme, ou 
une bombe camoufiée, est modélisée par une inclusion D C ^~ , de permitti- 
vité > et de perméablité /i^ > 0, enfouie dans la partie inférieure 
On supposera que D est suffisamment en profondeur par rapport à la taille 
de la couche, et que D est suffisamment grande devant les inhomogénéités du 
milieu inférieur. 

Notre problème modèle, certes simpliste, présente l'avantage d'être purement 
électromagnétique, et de prendre en compte la rugosité de l'interface. Il est 
nouveau ; pour le résoudre, nous proposons une démarche originale, en adap- 
tant pas à pas des techniques classiques. Son caractère multi-échelles et pério- 
dique conduit à l'approcher par un modèle effectif qui lui est équivalent, par 
des techniques asymptotiques usuelles p[T][5] . Par hypothèses sur l'objet D, 
son ajout dans le milieu périodique est équivalent à son ajout dans le miheu 
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effectif. La solution du problème initial en présence de D est alors calculée par 
la méthode des équations intégrales pi^ . dans le milieu effectif. 



1 Milieu effectif 



Le caractère multi-échelles périodique du problème ([T]) suggère d'appliquer les 
méthodes asymptotiques usuelles. La différence d'ordre des échelles permet 
d'effectuer l'étude en deux temps. Tout d'abord, de façon similaire à un calcul 
mené dans [2], on homogénéise le milieu inférieur ^~ en tronquant un dé- 
veloppement à deux échelles. Puis, de façon analogue à |T||5] . on élabore des 
conditions de transmission, de type impédance généralisée, permettant d'ap- 
procher l'effet de la couche mince. Celà s'obtient en écrivant un développement 
asymptotique à l'aide de correcteurs de couche limite décroissant exponentiel- 
lement. On trouve alors que la solution du problème modèle ([T]) peut être 
approchée par U solution d'un problème dans deux demi-espaces homogènes, 
avec des conditions de transmission à l'interface : 

■^AU + ùo^e-^U = 0, dans V • AVU + oo^e-U = 0, dans 



[U]=i^- Vt/|^+ , sur 70, 

^ d.,U\^-, - (0, 1)^ VUl- = V, ■ Vf/|^+' 



+^2 ■ {dl.,U\,^^ i-k^' - + ^3 U\^^ , sur 70. 



Les paramètres effectifs du milieu inférieur sont issus de l'homogénéisation, 
donnés par ^ = |^ Jy (I - Vxiv)) dy, et e' = {fj^ij + (1 - où 
X est la solution F-périodique à moyenne nulle du problème de cellule Vj; ■ 
fj_Y[y) ^ yX{y) = ^y' ij,Y{y) ^ ^ P*^^^ conditions de transmission, les coef- 
ficients sont donnés par : s = Ai, = ^î/^o, V^z = ^ - ^) Sq fiO, 1)"^, V^s = 

^u' (^^1 - e^^) lo' /, ^i/e = -sJr, ^olro + " ^) /r. ^olr, da, et 
(\i/o,V^o) est l'unique solution du problème de cellule : 

V ■ ^V^o = - ^) ^<^r, - - ^) ^-^ro, dans (0, 1) x M, 
î/i I — ^ ^o(yi, ■)> l-périodique, 
^'o — > 0, 2/2 -oo; *o — > ^0, 2/2 -> +oo, 

où par remise à l'échelle y := x/^, Ti et Tq sont les images respectives de 7ç 
et 7o, et À{y) = Als>-{x) + ^l^ci(x) + j^l^+{x). 

La solution U du problème effectif ([2]), et la fonction de Green G associée, 
peuvent être calculées par analyse de Fourier tangentielle à l'interface 70. Pour 
G, ceci nécessite alors l'évaluation d'intégrales de Sommerfeld. 
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2 Problême modèle d'objet dissimulé 



Nous avons transformé le problème modèle ([T]) en le problème effectif ([2]). 
Ajoutons une inclusion D C dans la partie inférieure du milieu, de per- 
mittivité 5/) > et de perméablité fiD > 0. On suppose que D est suffisam- 
ment en profondeur par rapport à la taille de la couche, et que D est grande 
devant les inhomogénéités du milieu inférieur. Sous ces hypothèses, l'analyse 
asjTiiptotique permet, formellement, de remplacer le milieu de fond par le mi- 
lieu effectif précédent. Pour résoudre ce problème, nous utilisons la méthode 
des équations intégrales. On note Gu la fonction de Green de l'anomalie ; 
on définit les potentiels de simple couche par S(p{x) := Jg^G{x,y)(p{y)da{y), 
Sd^{x) := /g^ Gd{x, y)ip{y)da{y). L'ajout de D transforme U enu par convo- 
lution avec les noyaux de Green [6|3fi] : 

u = U + Sip dans D^, u = S^il) dans 

où (fjip) G L?'{dD) X L^[dD) est solution d'un système d'équations intégrales 
permettant d'assurer les conditions de transmission sur dD. Ce dernier est 
résolu numériquement par la méthode des éléments finis de frontière. 
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